Chaos and plasticity in superconductor vortices: a low-dimensional dynamics 
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We present new results of numerical simulations for driven vortex lattices in presence of random 
disorder at zero temperature. We show that the plastic dynamics of vortices display dissipative 
chaos. Intermittency "routes to chaos" have been clearly identified below the differential resistance 
peak. The peak region is characterized by positive Lyapunov exponents characteristic of chaos, 
and low frequency broad-band noise. Furthermore we find a low fractal dimension of the strange 
attractor, which suggests that only a few dynamical variables are sufficient to model the complex 
plastic dynamics of vortices. 
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I. INTRODUCTION 

The physics of periodic structures in a random pinning 
potential has significantly improved these last years. In 
particular, an elastic approach of the static problem has 
brought conclusive results in many situationsi^^. 
When driven over the random medium, the fast mov- 
ing periodic structure is able to absorb dislocations that 
eventually appeared at intermediate velocities' 7 ^. In 
their elastic approach, Giamarchi and Le Doussal 9 show 
that a non-linear static disorder persists in the direction 
transverse to motion even at high velocity. They there- 
fore conclude that the moving lattice is a moving glass 
and not a moving crystal. As in the static case, in di- 
mension three and weak disorder the relative transverse 
displacements only grow logarithmically above a given 
length scale. The moving glass is therefore characterized 
by a quasi long range order with perfect topological or- 
der and rough static channels in which "particles" are 
flowing. In the case of point disorder, such moving struc- 
ture is called the moving Bragg glass. Dislocations be- 
tween the channels are expected for stronger disorder or 
in dimension two, creating decoupled channels and the re- 
sulting moving glass is called a moving transverse glass. 
It is characterized by a smecticlike order transverse to 
the direction of motion. However, the question of the 
stability of the fast moving phase is still controversial 10 . 
In the particular case of correlated disorder created for 
example by heavy ions irradiation in type II supercon- 
ductors, the vortex moving structure is a moving Bose 
glass^ characterized by a transverse critical force and a 
diverging tilt modulus due to the localization effect aris- 
ing from the columnar pins. This novel feature specific 
to correlated disorder results in the so-called dynamical 
transverse Meissner effect that has been confirmed in nu- 
merical simulations 12 . 

In the case of strong point disorder or intermediate 
velocities dislocations are likely to appear and plastic de- 
formations should be considered. In conventional super- 
conductors the vortex experiments show plasticity close 
to the melting line, i.e. in the peak effect region. In this 
region of the phase diagram, the shape of the current- 
voltage I — V curve is modified and displays a peak in 



the differential resistance dV/dl curv o 13 ^ 14 . Furthermore, 
the voltage noise measurements increase of several orders 
of magnitude in this regio n 14 : 15 . Recent studies explain 
the peak effect with the important role of surface pinning 
or surface barriers that lead to non-uniform current flow 
in the sample 1 ^. In numerical simulations, the plastic- 
ity is observed when the pinning is not too weak lead- 
ing to filamentary depinning made of plastic flow around 
pinned region o 7 ' 8 ' 17 . In this case the intermediate veloc- 
ity regime is accompanied by a change of curvature in the 
velocity-force curve, generating a peak in the derivative 
curve. This is similar to the peak of the differential re- 
sistance curve obtained in the peak effect region of type 
II superconductors. In such plastic phases topological 
defects proliferate in the vortex lattice, and a complex 
and apparently very disordered dynamics of vortices is 
growing in size. It is therefore clear that these regimes 
cannot be described by an elastic approach. The theoret- 
ical understanding of plastic depinning and plastic flows 
remains an open problem. Different coarse-grain mod- 
els (where the Larkin domains are the degrees of free- 
dom) are developed to describe plastic deformations: see 
ReffhSl and references therein. In particular, Marchetti et 
al^ show the existence of a tricritical point which sep- 
arates continuous depinning transitions at weak disorder 
from discontinuous depinning transitions with hystere- 
sis at strong disorder. These mean field results agree 
with the vortex simulations of Refl20l. However, numer- 
ical studies of a phase slip model 21 suggest the absence 
of hysteresis at the thermodynamic limit. In another nu- 
merical phase slip model, the plastic depinning for strong 
disorder appears to be a continuous transition^. 

The theoritical approach of plastic flows is difficult 
due to their intrinsic complexity. However, we recently 
showed 23 that a natural approach in terms of dissipa- 
tive chaos gives some new understandings of the vortex 
plasticity. In particular, a low dimensional dynamics has 
been evidenced for N v = 30 vortices. In the present pa- 
per we show the pertinence of the chaotic approach for 
larger systems, but more importantly we show that the 
low fractal dimension of the chaotic attractor remains un- 
changed despite a much higher dimensional phase space. 
The outline of this paper is as follows. In section |TT] we 
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describe our numerical model for line vortices in type II 
superconductors. In section IIIII we show in details the 
transition to chaos followed by the vortices in the plastic 
phase. Section IIVI characterizes the chaotic phase itself 
which is evidenced in a wide range of the driving force 
around the differential resistance peak. Positive Lya- 
punov exponents, broad-band noise at low frequency and 
fractal dimension of the strange attractor are analysed in 
details. The crucial conclusion of a low dimensional dy- 
namics is explained, and the influence of the number of 
degrees of freedom and dissipation are discussed in sec- 
tion LY1 



II. NUMERICAL MODEL 

We consider N v Abrikosov vortices driven over a ran- 
dom pinning background in the (x, y) plane. At T = 
the overdamped equation of motion of a vortex i in po- 
sition r, reads 



dri 



(i) 



where is the distance between vortices i 



and j, rip is 

the distance between the vortex i and the pinning site lo- 
cated at r p , and Vj is the 2D gradient operator acting in 
the (x, y) plane. The viscosity coefficient is rj, F L — F L ~k 
is the Lorentz driving force due to an applied current. 
The vortex-vortex pairwise repulsive interaction is given 
by a modified Bessel function 

U vv ( nj ) = 2e A v K (r l3 /X L ), 

and the attractive pinning potential is given by 



-a v e 



-{r ip /R p f 



In these expressions A v and a p are tunable param- 
eters, Xl is the magnetic penetration depth, and 
(q = (4>o/4ttXl) 2 is an energy per unit length. We 
consider periodic boundary conditions of (L x ,L y ) sizes 
in the (x, y) plane. All details about our method 
for computing the Bessel potential with periodic 
conditions can be found in ReflH. A Runge-Kutta 
algorithm for molecular dynamics simulation is used for 
N v = 30, 270 and 1080 vortices in a rectangular basic 
cell (L x ,L y ) — (5, 6v3/2)nAi, where n = 1, 3 and 6 
respectively. We choose the same density of pinning 
centers as the vortex density. We consider the London 
limit k = Xl,/£, = 90, where £ is the superconducting 
coherence length. £ appears in the inner cutoff which 
removes the logarithmic divergence of Ko(t/Xl) at r — > 
(see Ref. [2J). The average vortex distance ao is set to 
a = X L , and R p = 0.22 Xl, r] = l,A v = 2.83 x 10~ 3 A L . 
We computed different pinning strengths a p /A v ~ 0.05, 
a p /A v ~ 0.35 and a p /A v ~ 1.05, corresponding to a 
maximum pinning force of respectively F^ ax ~ 0.2F , 
F%x ~ 1AF o and F%> ax ~ 4F , where F = 2e ^„/A L 
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Figure 1: (Color online) Intermittency route to chaos char- 
acteristics observed for N v = 30 vortices and for a pinning 
strength a p /A v ~ 0.35 : a) part of the time evolution of the 
longitudinal velocity V^ m {t) obtained for F L = 1.116 x 10~ 3 . 
One sees laminar (i.e. periodic) phases interupted by chaotic 
bursts of large erratic velocity fluctuations, b) distribution of 
the laminar phase durations of V£ m (t). Both plots indicate 
type II intermittency (see text). 



is a force defined by the Bessel interaction. The driving 
force applied along a principal vortex lattice direction x 
is slowly varied from up to a value far from the critical 
Lorentz force along x. 

The successive regimes we observe are pinned configu- 
rations where all vortices have zero velocity, followed 
by plastic channels flowing through pinned regions. 
The successive depinnings of the pinned regions result 
in a complex plastic phase of interconnected flowing 
channels. Finally coupled or decoupled channels mostly 
aligned with the driving force appear in the high driving 
phase. 



III. INTERMITTENCY "ROUTE TO CHAOS" 

We first begin with the transition to chaos (often 
called routes to chaos in dynamical system theory) that 
we observe for vortices in a very short applied force range 
below the differential resistance peak. In particular we 
show the results we obtain for N v = 30 vortices and for 
the pinning strength a p /A v ~ 0.35. In Fig. [T]a we show 
the typical longitudinal velocity of the vortex center of 
mass V£ m (t) that we measure in time for a given Lorentz 
force F L = 1.116 x 10~ 3 - 0.3943F . V£ m (t) shows time 
intervals where the motion is almost periodic (laminar 
phases) and interrupted by chaotic bursts displaying 
large erratic velocity fluctuations. The almost periodic 
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motion corresponds to a plastic channel flowing through 
pinned regions, and the instability gives way to a chaotic 
burst where all vortices are moving erratically. Later 
the system goes back to the almost periodic regime with 
pinned vortices and moving vortices that synchronize 
temporarily their motion. Again the instability gives 
way to another chaotic burst, and so on. Therefore, such 
signal V£ m (t) shows intermittency which may be one of 
the three main known scenarios to drive a dissipative 
system from periodicity to chao o 25 ' 26 . In the case of 
intermittency, the dynamical system has a stable limit 
cycle below the intermittency threshold F t , whereas 
for F — ► the dynamical regime is intermittent with 
apparently periodic oscillations interrupted by large 
fluctuations. The amplitude and duration of these large 
fluctuations are almost the same from one fluctuation to 
the other, and they depend little on the force. When 
F — > F^~ these fluctuations become increasingly rare and 
disappear completely below the threshold. Therefore 
the mean frequency of the chaotic bursts goes to zero 
at the transition but not the amplitude neither their 
duration. The intermittency "route to chaos" has several 
characteristics and may be mainly classified in three 
types (I, II and III) depending on the way the limit 
cycle loses stability. In general a trajectory in phase 
space is linearly stable if all the eigenvalues (called 
the Floquet multipliers) of the Floquet matrix have 
magnitude smaller than 1, and loses stability if one of 
them leaves the unit circle in the complex plane 2 ^. A 
crossing by the (+1) value on the real axis generates a 
saddle-node bifurcation and gives type I intermittency. 
The crossing by (—1) on the real axis generates a a 
subcritical subharmonic bifurcation and gives type III 
intermittency, whereas type II intermittency is observed 
above a subcritical Hopf bifurcation and corresponds to 
a crossing of the Floquet multiplier by two conjugate 
complex values. To determine the type of intermittency 
corresponding to Fig. la, we first measure for a given 
value of the applied force the distribution of the laminar 
(i.e. periodic) phase durations. Fig. lb shows such dis- 
tribution obtained for V£ m (t) displayed in Fig. la. Such 
distribution is qualitatively very different from the one 
obtained for weaker pinning a p /A,, ~ 0.05 corresponding 
to F^ ax ~ 0.2fo shown in Refl23l, and characteristic 
of type I intermittency. The very different shape of 
the distribution we observe here for a p /A v ~ 0.35, in 
particular the long tail for large durations, cannot be 
attributed to type I intermittency but is expected for 
type II or type III intermittency. A possible way to 
discriminate between these two possibilities is to enlarge 
the signal V£ m (t) at the end of the almost periodic 
intervals, i.e. just before the chaotic bursts. Since type 
III intermittency is associated to a period doubling 
bifurcation (crossing by the (—1) value of the Floquet 
multiplier), the chaotic burst is expected to appear just 
after the increasing of a 1/2 subharmonic oscillation 
at the end of the laminar intervals. Such behavior is 
clearly not observed in the signal V£ m (t) displayed in 



Fig. la., and is confirmed by a power spectrum analysis. 
We therefore conclude that type II intermittency is the 
"route to chaos" followed by N v — 30 vortices in rather 
strong pinning a p /A v ~ 0.35 27 . Since type I intermit- 
tency was observed for weaker pinning a p /A v ~ 0.05 23 , 
we therefore show that the pinning strength may change 
the nature of the bifurcations brought into play. 



IV. CHAOS 

We now turn to the chaotic phase itself wich occurs in 
a large force range above the very narrow intermittency 
range. For the study of chaos we now concentrate on a 
larger system size of N v — 270 vortices flowing over a 
random medium made of N p — 270 pinning centers. We 
present results for a pinning strength of a p /A v ~ 1.05. 



A. Differential resistance 

Figure [5] displays the differential resistance curve Rd = 
dV£ m /dF L . The shape is very similar to the differen- 
tial resistance curve obtained in the peak effect region 
of type II superconductor a 13 i 14 . It shows a peak corre- 
sponding to the S shape of the velocity-force curve. In 
numerical simulations, this peak indicates plasticity in 
the vortex flowing. Above the main peak in Fig. [5] ap- 
pears the so-called fingerprint phenomenon already ob- 
served in experiments^^. This anomaly corresponds 
to the depinning of the last pinned vortices. Our study 
shows that for various system sizes and various pinning 
strengths, the plasticity of the filamentary vortex flow 
indicates chaotic dynamics. Chaos begins close to the 
bottom of the differential resistance peak where the in- 
termittency route to chaos is observed, and ends just 
above the fingerprint phenomenon. 



B. Lyapunov exponents 

An unambiguous signature of chaos is given by at least 
one positive Lyapunov exponent illustrating the sensi- 
tive dependence on initial conditions (SDIC), which is a 
property of chaotic attractors only. Starting from two 
neighbouring initial conditions on the chaotic attractor, 
the two corresponding trajectories in phase space will di- 
verge exponentially fast one from another while still stay- 
ing on the chaotic attractor. The Lyapunov exponents 
are the inverse of the characteristic times of this expo- 
nential divergence. Therefore, to compute the maximal 
Lyapunov exponent A one has to compute the distance 
d{t) between two initial neighbouring trajectories on the 
chaotic attractor. Since we integrate N v first order dif- 
ferential equations of motion (Eq. [T]) , the phase space is 
defined by the 2N V vortex coordinates and the distance 
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Figure 2: (Color online) Differential resistance Rd = 
dV£ m /dF L obtained for N v — 270 vortices and for a pin- 
ning strength of a p /A v ~ 1.05 . The main peak close to 
F L = 0.01 indicates the S shape of the velocity-force curve, 
and the so-called fingerprint phenomenon is observed between 
F L = 0.02 and F L = 0.03. As indicated by doted lines, chaos 
is measured between F L = 0.006 and F L = 0.03, as shown in 

Fig.rj 



d is defined by 



d 2 {t) 



N v 

E 

»=i 



Xi{t)-Xi{t)) +[Y i (t)-Y i (t) 



where X t (t) = x t (t) - x cr Jt), Y l (t) = y^t) - y cm (t), 
Xi(t) = Xi(t) - x cm (t), Yi(t) = yi(t) - y cm (t)- In 
these expressions, (xj, yi) and (xi, jji) are the vortex 
i coordinates, and (x cm , y cm ) and (x cm , y cm ) are the 
respective coordinates of the center of mass. The tilde 
notation (x, y) refers to the second trajectory generated 
by the neighbouring initial condition. Figure[3Ja displays 
two examples of the time evolution of d that we typically 
find in the region of the differential resistance peak. 
It clearly shows an exponential divergence d ~ e xt of 
the distance between two trajectories starting with an 
initial distance of 10 -12 between them (see the green 
curve). The slope therefore defines the positive maximal 
Lyapunov exponent A characteristic of chaotic dynamics. 
When the distance d becomes comparable to the size of 
the chaotic attractor, a saturation effect in the growth 
of d(t) naturally appears. Furthermore Fig [3]a shows 
that, for the same driving force, changing the initial 
distance between the two trajectories (see the orange 
curve) doesn't change the maximal Lyapunov exponent 
A (whithin the error bars), but also doesn't change the 
length scale d c below which the exponential divergence is 
observed, so that both A and d c appear as characteristic 
quantities of the intrinsic chaotic dynamics of the vortex 
system. Fig. [3Jb displays the evolution of A that we 



Figure 3: (Color online) Signature of chaos in the plastic 
phase for N v — 270 vortices and for a pinning strength of 
a p /A v ~ 1.05 : a) semi- log plot of the time evolution of the 
distance d(t) between two initial neighbouring trajectories in 
phase space for F L = 0.012. Two different initial distances 
between neighbouring trajectories are shown, 10~ 12 (green) 
and 10 -6 (orange). The positive slope A (see solid line as 
guide for the eye) indicates the exponential divergence of d(t) 
representative of chaos. The same slope observed up to the 
same length scale d c above which the growth of d saturates, 
indicate that both A and d c are characteristic quantities of 
the intrinsic chaotic dynamics, b) Evolution of the maximal 
Lyapunov exponent A (blue circles) with the Lorentz force, 
plotted together with the differential resistance curve Rd = 
dV^ m / dF L (red triangles). Each Lyapunov exponent is the 
average over 20 couples of initial conditions and the error bars 
are the standard deviation. 



measure in the chaotic region. In the same graph we 
plot the differential resistance curve Rd = dV£ m /dF L . 
The shape of the Lyapunov curve closely follows the 
differential resistance curve. In particular a peak 
appears in the Lyapunov curve simultaneously with 
the differential resistance peak, and a shoulder at 
F L ~ 0.2 indicates the fingerprint phenomenon shown 
in Fig. [2] The maximal Lyapunov exponent is in some 
sense a measure of the degree of chaos since the higher 
Lyapunov exponent A, the faster the divergence of two 
chaotic trajectories. Therefore Fig. [3]b shows that 
chaos is fully developed at the differential resistance 
peak. Furthermore, the shoulder in the Lyapunov curve 
close to F L — 0.020 indicates that the decreasing of 
chaos above the differential resistance peak is slowing 
simultaneously with the fingerprint phenomenon which 
confirms the idea of the depinning of the last pinned 
vortices. 
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Figure 4: (Color online) Power spectrum S x (f) of the longitu- 
dinal noise for N v — 270 vortices and for a pinning strength 
of a p /A v ~ 1.05: a) F L = 0.006, below the differential re- 
sistance peak. The lower dotted line corresponds to 1/rdiff 
while the higher corresponds to the maximal Lyapunov expo- 
nent A. b) F L — 0.009 corresponding to the peak of the differ- 
ential resistance. The dotted line corresponds to the maximal 
Lyapunov exponent A. In this case, Tdiff is larger than the 
experimental time window so that 1/rdiff does not appear in 
the power spectrum. The chaotic noises iV^ and N y shown 
in Fig. [5] are averaged in the range 1/rdiff < / < A. For 
F L — 0.006, note the colored noise signature S x (f) ~ f~"' of 
the diffusion process in the power spectrum for / < 1/rdiff. 



C. Low frequency noise 

Many experiments in type II superconductors have 
measured large excess noise, in particular in the peak ef- 
fect regio n 14 ' 15 or close to the first order transition where 
the vortex lattice melts in the vortex liquid phased, and 
very few show transverse noise measurements^. However 
the origin of the excess noise remains an open question 
since several mechanisms are proposed. In numerical sim- 
ulations, the differential resistance peak and excess noise 
originate from plastic flow of vortices. We analyse plas- 
ticity in terms of chaos and in particular we argue that 
the low frequency broad-band noise measured in the dif- 
ferential resistance peak region is a signature of the SDIC 
property which is characteristic of chaos since it implies 
the loss of memory for the chaotic system. 
We compute the power spectrum S a (f) of the velocity of 
the center of mass V£ m (t) defined by 



SM) 



i 



dt V£ m (t) exp{i2nft) 



where a — x or y. Figure [4] shows the typical power spec- 
tra we obtain close to the differential resistance peak. 
The broad-band noise at low frequency is obvious and 
shows the impossibility of long term prediction. How- 
ever, the washboard frequency / = V x m /ao appears in 
the power spectra immediately above the peak of the 
differential resistance (not shown) therefore showing the 
beginning of temporal order of the lattice although still 
chaotic. We shall now compute the broad-band low fre- 
quency noise characteristic of chaos. Two time scales ap- 
pear in the region of the differential resistance peak. 1 /A, 
where A is the Lyapunov exponent, is the characteristic 
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Figure 5: (Color online) N v = 270 vortices and for a pin- 
ning strength of a p /A v ~ 1.05 : differential resistance (green 
triangles), and longitudinal N x (red circles) and transverse 
N y (blue squares) low frequency noise averaged in the range 
1/Tdiff < / < A. 

time above which chaos appears and Tdiff is the char- 
acteristic time above which diffusive motions are mea- 
sured. We compute w x (t) = 7^X^=1 [Xi(t) — -^(O)] 2 
and w y (t) = [Yi(t) - YiiO)} 2 , with X t (t) = 

Xi(t)-x cm (t) andFi(i) = yi(t)-y cm (t). As already found 
by Kolton et al£ we find w x (t) ~ <^ x and w y (t) ~ t* v 
with exponents indicating normal and anomalous diffu- 
sion. We observe these diffusive motions for time scales 
larger than Tdiff, which varies with the driving force: a 
maximum of Tdiff is found and coincides with the dif- 
ferential resistance peak. Such diffusive motions have a 
clear signature in the power spectra S a {f). Indeed, for 
very low frequencies corresponding to time scales larger 
than Tdiff a colored noise S a (f) ~ / ' appears in the 
power spectrum (see figure 0^). Therefore, to measure 
the degree of chaos we compute the average noises N x 
and N y over the low frequency range 1/rdiff < / < A, 
i.e. 



1 



A - 1/Tdiff Jl/ 



S a (f) df 



T diff 



where a — x or y, and the Lyapunov exponent A have 
been evidenced in Fig. [3l Concomitantly with the differ- 
ential resistance, Fig. [5] shows the longitudinal N x and 
transverse N y low frequency noises averaged in this way 
in the chaotic region. Again the shape of the noise curves 
closely follows the differential resistance curve. The max- 
imun of the longitudinal noise coincides with the peak of 
the differential resistance which confirms that chaos is 
fully developed at the peak. Furthermore, anomalies in 
the low frequency noises appear simultaneously with the 
fingerprint phenomenon. Finally note that the maximal 
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longitudinal noise level is almost three orders of mag- 
nitude higher than the level measured at the end of the 
chaos region, i.e. just above the fingerprint phenomenon. 

D. Fractal dimension of the strange attractor, and 
low-dimensional dynamics 

In sections [IV B I and II V CI we characterized the absence 
of temporal correlation in the chaotic regime due to the 
SDIC fundamental property. We shall now character- 
ize the spatial correlations within the chaotic regime by 
computing the dimension of the chaotic attractor. The 
peculiarity of a chaotic attractor comes from its two 
properties that seem hard to reconcile. Trajectories of 
phase space converge and remain confined to a bounded 
region (attraction), although neighbouring trajectories 
on the attractor separate exponentially fast (SDIC). It 
involves a double mechanism: stretching necessary for 
the SDIC property, followed by folding required for con- 
fined trajectories to a bounded region. Repeated applica- 
tions of stretching and folding generate a fractal attractor 
whose peculiar properties justify the name strange attrac- 
tor. There exist several definitions of a fractal dimension 
which are not equivalent and lead to different numeri- 
cal values for a given object. The Hausdorff-Besicovitch 
dimension is probably the most famous but rather un- 
workable in practice. A more efficient method is to com- 
pute the correlation dimension proposed by Grassberger 
and Procaccia^i. Consider a set of very many points 
{Xi(t), Yi(t)} (see section HVBj) on the strange attractor 
generated by letting the system evolve a long time. The 
correlation dimension v is the exponent of the power law 
C(p) ~ p v , where 

^ m 

C(p) = linim^oo — ^ H (P - Pkl) 

m k,l=l 

measures the number of couples of points (fc, I) on the 
chaotic attractor whose distance pu is less than p. H(z) 
is the Heaviside function. To estimate the correlation 
dimension v, we determine the local slope on a log-log 
plot of C{p) defined by 

_ d{log 10 C(p)) 
d{log w p)) 

v is deduced from the intervals where v\ oc is constant, i.e. 
where a true power law is measured. Figure [6] displays 
vice with respect to logiop that we obtain for F L — 0.008. 
It clearly shows a constant value of vi oc ~ 1.12 over a lim- 
ited range of p. It shows the validity of the power law 
behavior C(p) ~ p v in this range. The correlation di- 
mension for F L = 0.008 is therefore v ~ 1.12. Above 
and below the limited range where the power law behav- 
ior is observed, two regimes exist that are well known 
for such correlation dimension computations. Because 
of the finite number of points {Xi(i), Yi(t)} defining the 
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Figure 6: (Color online) Local slope vi oc with respect to 
logio(p) for F L = 0.008 in the case of N v = 270 vortices and 
a pinning strength of ot v /A v ~ 1.05 . A constant value of ui oc 
appears showing the validity of the power law C (p) ~ p" . The 
value v ~ 1.12 for the correlation dimension is measured for 
F L — 0.008 in the power law regime. Below and above this 
regime, are two regimes that clearly deviate from the power 
law (see the text for details). 

strange attractor in our simulations, the number of pairs 
(k,l) whose distance is less than p is small for small val- 
ues of p, giving therefore a poor statistics. And above the 
limited range where vi oc is constant, p becomes compa- 
rable to the size of the strange attractor, and increasing 
p will not add new (k, I) couples. Therefore C(p) — > 1 
and such saturation leads to vi oc — > 0. 
The computation of the correlation dimension in the 
whole chaotic region, i.e. in the region of the differential 
resistance peak, shows the important result that 

1 < v < 1.5 . 

This result indicates that in our simulations where the 
time is discrete the chaotic trajectories in phase space re- 
main confined in a subspace of dimension less than two. 
Since the discrete time system can be considered as a 
Poincare section (dimension D—l) of the continuous time 
system (dimension D), we conclude that the continuous 
time system of differential equations (Eq. [T|) has chaotic 
trajectories confined in a subspace of dimension less than 
three. This is a crucial result since it indicates that the 
complex plastic vortex dynamics of very many degrees of 
freedom is low dimensional, and can be modeled by only 
three dynamical variables. 

Note that the fractal dimension of the chaotic attractor 
that we find is strikingly similar to the fractal dimen- 
sion of the vortex channel network in the 2D real space 
found in Ref. [32| for different pinning strengths. Such 
coincidence is not clear since the link between chaotic 
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trajectories in phase space and in real space is inaccessi- 
ble. 



V. DISCUSSION 

The same result 1 < v < 1.5 was obtained for much less 
degrees of freedom (60 instead of 540 considered here) in 
our previous paper—. We also checked that for the much 
larger system N v = 1080, i.e. 2160 degrees of freedom, 
the chaotic trajectories remain bounded in a low dimen- 
sional fractal attractor with 1 < v < 1.5. Therefore the 
important result of a low dimensional dynamics of the 
continuous time system of differential equations Eq. [T] on 
a chaotic attractor of dimension less than three remains 
valid whatever the number of degrees of freedom Ndf is. 
Since the system size does not seem to play any role in 
the low dimension of the strange attractor, we investigate 
the role of dissipation which shrinks the volumes in phase 
space and might influence the dimension of the attractor. 
Indeed, the dissipation of a dynamical system is given by 
the divergence of the flow. Since Eq. [T] is a first order 
differential equation of the form f = F(r) the divergence 
V.F of the flow is not controlled by the viscosity coef- 
ficient r] but by the vortex-pin and vortex-vortex force 
component derivatives ^ i {dF i /dx i + dFi/dyi). Keep- 
ing the same parameters £, \l and A v for the line vor- 
tices, we can therefore modify the dissipation of the flow 
by modifying the pinning parameters, e.g. the pinning 
strength a p . We therefore computed the fractal dimen- 
sion of the chaotic attractor for a p /A v ~ 0.05, i.e. for 
a pinning strength 20 times less. In this case the pin- 
ning is low since while increasing the driving force for 
N v = 30 vortices we observe a clear transition between 
chaos (plasticity) and the coupled channels regime that 
might be called the moving Bragg glass. Our compu- 
tations show that in the chaotic phase, the same result 
1 < v < 1.5 appears. Furthermore, for the same low 
dissipative flow where a p /A v ~ 0.05, and for a much 
larger system size N v = 1080 vortices, the same result 
1 < v < 1.5 holds. These computations therefore clearly 
show that the low fractal dimension of the chaotic attrac- 
tor that we measure remains unchanged for larger system 
sizes 30 < N v < 1080 or/and for lower flow dissipation. 



A natural extension of our present study is to anal- 
yse the chaotic dynamics of driven vortices in 3D. Pre- 
vious 3D studies 33 obtained different dynamical regimes 
depending on the coupling strength between planes and 
the driving force. In particular, it would be very inter- 
esting to compute the fractal dimension of the chaotic 
attractor in the 3D plastic and 3D smectic phases. 



VI. CONCLUSIONS 

In this paper we investigated the properties of plastic 
flows of superconductor vortices driven over a random 
medium. We interpreted the plastic dynamics in terms of 
dissipative chaos. The link with the differential resistance 
peak is established. Chaos begins at the beginning of the 
peak and ends after the "fingerprint" phenomena that 
occurs above the peak. We first identified with details 
the transition to chaos ("route to chaos"). The scenario 
we have measured in our numerical simulation is the in- 
termittency. Type I and type II intermittency occur de- 
pending on the pinning strength. We then studied in de- 
tails the chaotic dynamics that occurs in the force range 
of the differential resistance peak. Positive Lyapunov ex- 
ponents are measured which is a non-ambiguous signa- 
ture of chaos. A low frequency broad-band noise has also 
been evidenced in the chaotic region. Both Lyapunov 
exponents and low frequency broad-band noise closely 
follow the evolution of the differential resistance curve, 
in particular the existence of a peak. Finally, we mea- 
sured the fractal dimension of the chaotic attractor. We 
showed that the discrete time system in our simulation 
is confined to a subspace of dimension less than 2. We 
therefore conclude that the continuous time system of 
differential equations Eq. Q] is confined to a subspace of 
dimension less than 2 + 1 = 3. It suggests that the plas- 
tic dynamics of vortices may be described by a model 
with only three dynamical variables. We showed that 
this important result still holds for systems as large as 
2160 degrees of freedom and also for systems with low 
dissipation. Our results give new understandings of the 
plastic flow phase and open new perspectives for further 
theoritical studies of plasticity. 
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